Accepted for publication in ApJ Letters 

Preprint typeset using I^T^X style cmulatcapj v. 3/16/03 



o 
o 

(N 
Or 

(N 



STRONG GRAVITATIONAL LENS CANDIDATES IN THE GOODS ACS FIELDS 1 

CD. Fassnacht 

Department of Physics, University of California, 1 Shields Avenue, Davis, CA 95616 

AND 

L. A. Moustakas, S. Casertano, H. C. Ferguson, R. A. Lucas, Y. Park 2 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218 
Accepted for publication in ApJ Letters 

ABSTRACT 

We present results from a systematic search for strong gravitational lenses in the GOODS ACS data. 
The search technique involves creating a sample of likely lensing galaxies, which we define as massive 
early-type galaxies in a redshift range 0.3 < z < 1.3. The target galaxies are selected by color and 
magnitude, giving a sample of 1092 galaxies. For each galaxy in the sample, we subtract a smooth 
description of the galaxy light from the zsso-band data. The residuals are examined, along with true- 
color images created from the -6435^06*775 data, for morphologies indicative of strong lensing. We 
present our six most promising lens candidates, as well as our full list of candidates. 
Subject headings: galaxies: general — galaxies: high-redshift — gravitational lensing — surveys 
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1. INTRODUCTION 

A great strength of gravitational lenses as astrophysical 
tools is that they provide direct measures of total mass, 
without the requirement that the mass be luminous or 
baryonic. At the most basic level, the image separa- 
tion in lens systems provides a nearly model-independent 
measuremen t of the projected mass of the lensing 
galax y (e.g., lKochaneklll991t iWambsganss fc Paczvnskil 
1994). The combination of additional observational 
constraints and more sophisticated modeling can pro- 
vide data on the elongat ion and orientation of the 
mass distribution (e.g., iKeeton. Kochanek. fc Seliakl 
Il997|) . the radial mas s profile of the lensing galaxy 
(e.g. , IKochanek et al.l I200TI IKoonmans fc TreiJIF 



IRusin et al.l2002UWucknitz et al 120031) . an d the possible 
existence of dark matter subhaloes (e.g., |Brjyia5 eLaL| 
1 20021 iDalal fc Kochanekl 120021: iMao fc Schneider! 119981 
iMetcalf fc Madaul !200*lT iMoustakas fc MetcaU 120031 . 
Therefore, gravitational lenses provide a wealth of in- 
formation about moderate redshift galaxies and can 
provide important insights into ga l axy evolution (e.g . , 
IKeeton. Kochanek. fc Seliakl 11991 IRusin et alJ [2003). 
An increase in the number of known lenses benefits al- 
most all these studies and thus motivates searches for 
new lenses. The data obtained with the Advanced Cam- 
era for Surveys (ACS) on the Hubble Space Telescope 
as part of the Great Observator ies Origins Deep Survey 
(GOODS; Giavalisco et al. 2003) offers an unprecedented 
combination of angular resolution, depth, and sky cover- 
age. In this paper we report the results of a systematic 
search for lenses in the data obtained in the first three 
epochs of observations of the GOODS fields. 
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2. THE TARGET POPULATION 

The search strategy that we used was focused on 
the likely lens galaxies rather than the likely back- 
ground sources. This is in contrast to the searches con- 
ducted at radio wavelengths (e.g.. iBrowne et alJ l2003t 
iMvers et all2003l: iWinn et al.l2000|) and optical sea rches 
targeting known quasars (e.g.. iBahcall et al.lll992T) . For 
this search, we selected the objects most likely to pro- 
duce strong lenses, namely massive early-type galaxies 
at 0.3 < z < 1.3. Theoretical studies indicate that these 
galaxies should dominate the lensing population, with 
only ~10% or fe wer of strong lenses produced by spi- 
ral galaxies fe.g.. iFukugita fc Turnerlll991t iMaoz fc Ri^d 



H993HTurner. Ostriker. fc Gottlll984() . Similarly, a ray- 
tracing analysis con ducted on the norther n Hubble Deep 
Field data CHDF-N: IWilriams et allll996|1 indicated that 
nearly all the lensing cross-section was provided by mas- 
sive ellipticals at z ~ 1 (Blandford, Surpi, & Kundic 
2001). Thus, our target sample was chosen to maximize 
the chance of finding massive early-type galaxies in the 
desired redshift range. 

All candidate selection was done using the ACS data 
obtained in the F435W, F606W, F775W, and F850LP fil- 
ters (hereafter, B435, V606, *775j and zs5o)- We consider 
only objects with isophotal magnitudes Z850, iso < 23.0. 
Colors were measured in matched apertures, with the 
apertures defined by the Z85o-band data 3 . In a color- 
color plot of (S435 — 1775) vs (1775 — ^85o)j galaxies at 
z > 0.3 with old stellar populations should lie in the 
upper part of the diagram (see Figure 1). We define 
our target region as (-B435 — 1775) > 2.0, which pro- 
duces a sample of 1092 galaxies. Figure 1 also contains, 
as examples, no-evolution tracks for different types of 
galaxies, which show solely the effect of spectral features 
moving into different filter bands for galaxies at differ- 
ent redshifts. Although the track for elliptical galaxies 
in Figure 1 would suggest that our target region con- 
tains ellipticals in a wider redshift range than desired, 

3 For complete details of the source ext r action , photometry, and 
limiting magnitudes, see IGiavalisco et al.l <2003T) . 
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Fig. 1. — Color-color diagram for galaxies in the GOODS south- 
ern field 2g5o-selected catalog, limited to 2850, iso < 23.0. The tar- 
get galaxies for the lens search lie above the dashed line. The open 
symbols connected by lines represent the expected change with 
redshift of colors for non-evolving elliptical (E), Sb, and starburst 
(SB6) galaxies. The points are separated by Az = 0.1, starting at 
2 = and with solid points every Az = 0.5. On each track, the 
points corresponding to z = 0.5 and z = 1.0 are marked. 



the combination of the low volume probed at low red- 
shifts, cosmological dimming at high redshifts, and evo- 
lution implies that the redshift range selected is closer 
to the desired 0.3 < z < 1.3. Based on comparisons 
with both non-evol ving and passive evolution models 
l)Kinnev et al]ll99fih . as well as spectroscopy of ~ 200 
galaxies in the southern GOODS field, we find that the 
regions marked in the figure do contain the vast majority 
of the massive early-type galaxies in the relevant redshift 
interval. The target region also includes a large contam- 
ination (~ 20%) from highly-obscured systems (partic- 
ularly dusty bulge-dominated galaxies). Since our goal 
is to be as inclusive as possible, without needing to ex- 
amine the full catalog in detail, this gives an adequate 
sample of target galaxies. 

3. LENS- SEARCH TECHNIQUES 

3.1. Subtraction of galaxy models 

Our primary technique for searching the target sample 
for gravitational lens candidates was to model the light 
distribution of the target galaxy and then subtract the 
model from the data. The residuals for each target were 
examined for indications of gravitational lensing. To in- 
crease the robustness of our results, we used two methods 
to subtract the galaxy emission for each target. 

The first method used customized IDL scripts that pro- 
vided an empirical description of the target galaxy emis- 
sion. The galaxy centroid, ellipticity, and position angle 
were calculated from the light distribution. The input 
image was split into a series of concentric ellipses, each 
with the input ellipticity and position angle. The em- 
pirical description of the galaxy emission was then pro- 
duced by boxcar-smoothing the distribution of flux as a 
function of semi-major axis. This simple approach pro- 
vided acceptable representations of most of the target 



galaxies, and non-symmetric features were clearly visi- 
ble in the residual images produced by subtracting the 
smoothed image from the input. To judge the effective- 
ness of the technique in subtracting smoothly-distributed 
galaxy light, we use as a figure of merit the ratio of the 
counts in the residual image to the counts in the origi- 
nal image, summed within a 3"-diameter aperture. For 
^75% of the targets, this ratio was less than 0.05, in- 
dicating a good fit to the galaxy emission. All targets 
producing higher ratios had either close neighbors or the 
kind of asymmetric structure that the routine was de- 
signed to find, e.g., knots of star formation, bars, dust 
lanes, possible lensed images, etc. Thus, this method is 
sensitive to compact lensed images and to tangentially 
stretched arcs. However, complete or nearly complete 
Einstein rings would be subtracted away by this method. 
Therefore, we used a second galaxy-subtraction method 
to complement the empirical technique. 

The second technique involved fitting parametric mod- 
els to all candidate lens galaxies. Experimentation with 
various combinat ions of parametric laws (using GALFIT; 
iPeng et"aT]l2002() showed that a single Sersic fit did re- 
markably well for the majority of galaxies in our sample, 
with a median Xfu/ 1 ' ~ 0.9. In most of the cases, the 
residuals produced by the GALFIT subtraction did not 
differ qualitatively from those produced by the empirical 
method, and no believable lens candidates were found 
solely from an examination of the GALFIT residuals. 

3.2. Selection Criteria 

The residuals produced by the galaxy-subtraction 
methods were carefully examined, as were the unsub- 
tracted true-color (-B43sV(306«775) images of the target 
galaxies. The latter provided invaluable color informa- 
tion about features in the target images, which was incor- 
porated into the candidate selection process. The crite- 
ria used to select lens candidates were based on features 
commonly seen in HST images of strong gravitational 
lenses 4 . Any one of the following features seen in either 
the residual images or in the true-color images would 
cause a target to be flagged as a possible lens candidate: 
(1) two compact images on opposite sides of the target 
galaxy; (2) multiple (more than two) compact images ar- 
ranged in a typical lensing morphology, (3) an arc, pos- 
sibly with a faint counterimage, tangentially stretched 
relative to the target galaxy, especially if the arc is bluer 
than the galaxy, (4) a ring of possibly blue emission sur- 
rounding the target galaxy. Because gravitational lens- 
ing is an achromatic phenomenon, the alleged lensed im- 
ages should have the same colors. However, some leeway 
was allowed to incorp orate the possible effects of differ- 
entia l extinction (e.g.. lBlandford et alJl200ll:lFalco et alJ 
1999). Two of the authors independently examined all 
of the target galaxies. The list presented in Tabic ^ is 
the union of the two selections. 

4. DISCUSSION 

The systematic search produced a list of 48 candidates 
of varying quality. The biggest difficulty in selecting can- 
didates was to differentiate rings or arcs of star forma- 

4 The CfA-Arizona Space Telescope Lens Survey web site, 
at http://cfa-www.harvard.edu/castles/ provides an excellent 
compilation ot such images. 
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Fig. 2. — Plots showing the B435 V606*775 images for our top six lens candidates. All plots are 3" on a side with the exception of the 
plots for J033213-274946, which is 4.5" on a side, and J123730+621301, which is 7" on a side.. 



tion from lensing morphologies. Most of the candidates 
have morphologies consisting of arc- or point-like blue 
features near red galaxies, with no clear counterimages; 
these were accepted as candidates because the faint coun- 
terimages might be below the detection threshold. All of 
the authors voted on which of the 48 candidates were 
most likely to be lenses. We list our six top candidates 
in Table [21 and show them in Figure [21 Further fol- 
lowup is necessary to establish more definitively which 
of the systems are real gravitational lenses. Therefore, 
we present the full list of candidates and their coordi- 
nates in Table ^ so that interested readers can make 
their own selection of most promising lens candidates 
and target selected systems for additional observations 5 . 
We note that the northern GOODS field includes the 
HDF-N area, which was searched extensively for gravi- 
tational lenses. The two best lens candidates found in 
those searches, J123652+621227 and J123656+62 1221 
ijHogg et all 119961: iZepf. Moustakas. h David Il997|) . are 
detected in the GOODS data. The former is fainter than 
our magnitude cutoff at Zs5o,iso — 23.7; the latter is in 
the target sample and was independently selected as lens 
candidate GDS J123657+621220 (Figure [2J. 

In this letter, we will provide only first-order estimates 
of the properties of our six best candidates. Systems for 
which later observations clearly establish the lensing na- 
ture will be treated more fully in forthcoming papers. A 
nearly model-independent quantity that may be derived 
from lens system data is the mass of the lensing galaxy 

5 See also http : //www. stsci . edu/science/goods/lensing/Lensc 

for images ol the full list oi candidates. 



within a cylinder of radius 9e, the Einstein ring radius. 
In cases with clear multiple images, a good estimate of 
9e is simply half the image separation. However, in most 
of the lens candidates presented in this letter, an arc is 
seen, but no counterimage is detected. Although this 
presents a difficulty, we note that in cases where an arc 
is highly magnified compared to the counterimage, the 
arc lies very close to the critical curve of the lens. Thus, 
in these cases the distance from the galaxy to the arc 
provides a reasonable estimator of 9e- We have used the 
F-band images to estimate Oe for the candidates because 
both the blue and red objects are detected at high signif- 
icance in these images. In Table [21 we provide estimated 
Einstein ring radii (0E,est) as well as the target galaxy 
magnitudes in apertures with radii 9 — 0E,est- More ex- 
tensive descriptions of the properties of the lens galaxies 
will become possible with spectroscopic data and more 
detailed modeling. 

5. SUMMARY AND FUTURE WORK 

In this paper we have presented the results of a system- 
atic search for strong gravitational lens candidates in the 
GOODS ACS data. The search targeted massive early- 
type galaxies in the redshift range 0.3 < z < 1.3, i.e., the 
sources most likely to produce strong lenses. We have 
presented a list of 48 lens candidates as well as images 
of our six most promising candidates. The search was 
conducted on images that are not yet at the full depth 
expected from the GOODS survey; however, we do not 
expect the ^30% increase in sensitivity provided by the 
ands/toaaEpfiira.1iieBiiL 3f the last two epochs of the GOODS ob- 
servations to yield a significant increase in the number of 
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Table 1. Lens Candidates 



ID 




RA(2000) 


Dcc(2000) 


B435i3o 


V606, iso 


^775, iso 


2850, iso 


Zphot 


CDS J033206- 


-274729 


03 32 06.431 


-27 47 28.76 


23.83 


22.94 


21.78 


21.09 


0.96 


GDS J033211- 


-274650 


03 32 11.403 


-27 46 49.98 


24.42 


22.77 


21.70 


21.34 


0.46 


GDS J033213- 


-274946 


03 32 13.006 


-27 49 46.08 


24.02 


21.98 


20.65 


20.19 


0.55 


GDS J033215- 


-274157 


03 32 14.825 


-27 41 57.17 


25.38 


23.60 


22.37 


22.00 


0.55 


GDS J033216- 


-274714 


03 32 15.805 


-27 47 13.61 


24.44 


22.44 


20.98 


20.51 


0.62 



Note. — Units of right ascension are hours, minutes, and seconds. Units of declination are degrees, arcminutes, and arcseconds. [The 
complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample.] 

a Photometric redshifts determined as described in|Mobashcr ct al. (2003). Photometric redshifts for the northern GOODS sources are 
being determined. 



lens candidates. The best candidates are being targeted 
for spectroscopic observations, in order to confirm or re- 
ject the lensing hypothesis. Based on these data, we will 
be able to measure accurate masses, mass-to-light ratios, 
and other parameters of the confirmed lens systems. 
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Table 2. Parameters of Candidate Lens Galaxies 



ID 




RA(2000) 


Dec(2000) 


Zphot 


d E ,est h 


B 4 S5,E C 


Vk>6,s c 


*775,B C 


2850, S C 


Notes d 


GDS J033213- 


-274946 


03 32 13.006 


-27 49 46.08 


0.55 


0.96 


24.82 


22.35 


20.93 


20.48 


Blob + possible CI 


GDS J033231- 


-275612 


03 32 31.067 


-27 56 12.34 


0.57 


0.62 


25.00 


23.10 


21.82 


21.42 


Arc 


GDS J033240- 


-274910 


03 32 39.598 


-27 49 09.57 


0.75 


1.31 


25.52 


23.18 


21.69 


20.80 


Faint CI 


GDS J1236414 


-621125 


12 36 41.139 


+62 11 25.32 




0.73 


26.82 


24.95 


23.57 


22.73 


Faint arc 


GDS J1236574 


-621220 


12 36 56.642 


+62 12 20.45 




1.00 


25.77 


24.25 


22.71 


21.97 


Arc 


GDS J1237304 


-621301 


12 37 29.900 


+62 13 01.15 




1.98 


23.10 


21.29 


20.22 


19.85 


Multiple arcs (+CI?) 



Note. — Units of right ascension are hours, minutes, and seconds. Units of declination are degrees, arcminutcs, and arcseconds. 
a Photometric redshifts determined as described in Mobashcr ct al. (2003). Photometric redshifts for the northern GOODS sources are 
being determined. 
b Estimated Einstein ring radius, in arcseconds. 
c Magnitude within an aperture of radius 6 E est- 
d CI = counter-image. 



